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The study of capillary wave scattering by a circular region with diflerent interfacial properties 
from the rest of an otherwise homogeneous interface is motivated by experiments on wave attenuation 
at a monolayer covered air-water interface where domains of one surface phase are dispersed in a 
second surface phase. Here the scattering function is calculated for an incident wave of frequency 
1/^ ' ll> (wavevector ko) scattering from an isolated circular domain of radius a with surface tension ai 

^\ , which is imbedded in an otherwise infinite interface of surface tension ctq. The underlying fluid is 

0^ ' treated as irrotational and the three-dimensional flow problem coupling the heterogeneous surface 

' to the underlying liquid is reduced to a set of dual integral equations, which are solved numerically. 

f— I , With this solution the scattering amplitudes and the total scattering cross sections are calculated as 

' a function of the surface tension ratio ao/(Ti and incident wavenumber koa. The analogous problem 

' of a discontinuous change in bending rigidity is also considered and the solution to the complete 

. viscous problem is outlined in an appendix. Experimental implications of these results are discussed. 

^ . I. INTRODUCTION 

> ■ 

' The interfacial region separating two distinct fluids has significant dynamical effects on many transport processes. 
, In particular, fluctuations in the interface shape are present in problems as diverse as classical propagation of capillary- 
^-H ■ gravity waves, fluctuations of cell membranes, and mass transport across the ocean surface. The majority of analyses of 
] such surface waves have treated a fluctuating interface with homogeneous properties. Nevertheless, due to the typical 
phase-separate|d|-structure of monolayer covered interfaces or the presence of impurities, the surface may in fact be 
heterogeneous .liH For example, interfacial tension or bending rigidity, two material properties used to characterize 
, the interfacial resistance to deformation, may vary along an interface. In this paper, we analyze the effect of a 

a circular domain of different interfacial properties on the propagation of a capillary wave and so deduce the scattering 
I , characteristics. 

'"^ ' A typical feature of surface fluctuations with wavelength A is that the velocity in the surrounding fluid may be 
' disturbed to a distance A with a magnitude proportional to the magnitude of the restoring force. Hence, in the 
O neighborhood of a substantial change in surface properties, it is likely that local transport processes are affected. A 
^ first approach to understanding the dynamics of a heterpgeneous interface is to describe the mechanism of a plane 
^ ■ wave scattering from an isolated surface inhomogeneity.u However, almost all previous experimental and theoijetical 
■ research on surface wave scattering has focused on the effects-of solid bodies in contact with the bulk fluidp e.g. 

scattering producetLby nonuniform bottomsOij, floating objectaZI'u, and solid structures such as cylinders protruding 
' from the interfaceE^EI These studies have concentrated on the gravity wave limit with a few also treating capillary wave 
propagation. There is a rich mathematical literature describing gravity wave scattering by floating surface objects 
(such as a circular disk) with papers dating back 40 years.El 

There are now many experimental observations, including surfactant systems and lipid monolayers, of the co, 
plex, heterogeneous surface structures that develop when amphiphilic molecules reside at an air-water interface 
Frequently, the interface consists of circular domains (with radii tens of microns ranging up to a centimeter) of one 
phase interspersed in a second phase with different surface concentration. The surfactant concentration differences 
locally change the parameters which govern the stress balance at the interface. Both the surface tension and bending 
rigidity of the interface can thus be altered and in these small mesoscopic systems, it is surface tension, or bending 
rigidity, rather than gravity, that provides the dominant restoring force to surface deformations. Therefore, we study 
capillary wave propagation in the presence of an isolated circular island with different surface properties from the 
surrounding interface. Our work provides an example of wave scattering due to a boundary condition heterogeneity. 

Experiments on capillary wave propagation at fluid-gas interfajces covered with surfactants can show a dramatic 
increase in the damping rate of the surface capillary waves. l2I"E3 In these experiments, a wave maker generates a 
capillary wave and either the dispersion relation, wave amplitude, or the width of the spectral peaks generated 
from light scattering of a thermally fluctuating interface, is measured, and the damping coefficient characterizing the 
attenuation of the interface deformation is extracted. The damping coefficient is observed to depend on the surface 
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coverage and appears to have a maximum near the regime where surfacejphases coexist. t3~ll3 The damping rate may 
also depend on the typical domain sizes which coexist in the monolayer .ciJ Since the damping rates may be five times 
greater than that due_to viscous damping, a mechanism including the effects of wave scattering by surface domains 
has been suggested,aE3 but there is as yet no proven explanation for the substantial damping rate increase. 

Boundary conditions for gravity waves do not involve higher mixed derivatives and as gravity acts uniformly along 
the interface, surface wave scattering may be produced by imposed boundary conditions, e.g. a fixed dock occupying a 
half-plane on an infinite free surface. Such problems have been studied with Wiener-Hopf techniquesO Furthermore, 
wave scattering from an intrinsic boundary condition heterogeneity can occur when a varying interfacial tension 
or rigidity is present. In this case, the higher mixed derivatives that occur in the boundary conditions couple to 
the nonuniform surface parameters to produce surface wave scattering. This type of problem has received little 
attention. In particular, approaches to this latter problem in the past have been approximate,a-pr treated a specific 
one dimensional geometry with additional assumptions concerning the stability of the problem.EJ 

Here, we study analytically surface wave scattering from a circular discontinuity in surface tension or bending 
rigidity. In Section II, we outline the boundary value problem for an irrotational subphase. Hankel transforms may 
be used to describe the velocity potential and application of the normal stress balance at the heterogeneous interface 
leads to a system of dual integral equations with Bessel function kernels. This system is solved numerically. The 
analogous problem accounting for viscous effects is described in an appendix. For a benchmark comparison, we also 
derive results for scattering from a stretched, free membrane with an imbedded region of different tension. In Section 
III, we consider capillary wave scattering from a circular region of different interfacial bending rigidity and compare the 
results to bending wave scattering along a membrane with a circular inclusion of different flexural rigidity. The results 
of these different models are discussed and compared in Section IV. A discussion of the experimental applicability of 
these two problems, the criteria for neglecting viscous effects, and the implications for capillary wave damping are 
also given. 



II. CAPILLARY WAVES AND SURFACE TENSION DISCONTINUITIES 



A. Formulation 



Consider the propagation of small amplitude capillary waves on a fiuid with kinematic viscosity v, density p and 
of depth H (Figure la). A plane wave from far away impinges on a circular domain of radius a which resides on 
the otherwise homogeneous interface. First, we treat the classical capillary wave problem with surface tension as 
the restoring force. Second, we study scattering from a bending rigidity discontinuity in Section III. In all cases, 
gravitational effects on wave propagation are neglected. 

For a surface tension discontinuity, the constant surface tension of the interface exterior to the domain is denoted 
tJo and of the domain is denoted ui. Mechanical stability of such a circular domain requires that a line tension act 
along the domain bouudary. In monomolecular films, the line tension arises from the van der Waals interactions 
between the molecules. cla For the dynamical calculation reported here, only effects of surface tension are important 
and line tension effects are higher order since the domain is assumed to deviate only slightly from a circular shape. We 
thus apply linearized capillary wave theory to obtain the scattering produced by the discontinuous region of surface 
tension. Both the velocity field and height of the fluctuating iatpjf ace are calculated. This interface displacement is 
often measured experimentally by light scattering techniques.li3iij 

Assuming an irrotational flow, the velocity may be represented by a potential 



^.(..t).^^. (2.1) 



The flow is assumed to be incompressible so that 

V20 = O. (2.2) 

It is convenient to use cylindrical polar coordinates (i?, 6, Z) and later the dimensionless coordinates (r, 9, z) are used. 
For small amplitude disturbances, the momentum equation is linearized, so that the pressure p = —pdt4>, and the 
boundary conditions are also linearized and applied at Z = 0. Finally, the surface height is denoted ri{R,6,t). Effects 
of viscosity are treated in Appendix A. 

The boundary conditions to be applied are a kinematic constraint 
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a dynamic condition which follows from the normal stress balance, 



dt 



at Z = 0, 



(2.4) 



where V^^ is the two-dimensional Laplacian expressed in the surface coordinates, and an impenetrability condition at 
the bottom of the trough 



at Z = -H. 



Taking a time derivative of (2.4) and using the kinematic condition (2.3) yields 



(2.5) 



(2.6) 



which, along with (2.2) and (2.5), govern the evolution of the velocity potential. 



For a plane wave propagating in a fluid with homogeneous surface properties and depth H, 



^ikoX — iujt 



H) determines the well-known dispersion relation for surface capillary waves excited at frequency w. 



coshfe(Z- 



puj 



crfcg tanh kfjH, 



(2.7) 



where fcg = 27r/A denotes the wavevector and A is the wavelength. The normal stress contribution to (2^) due to the 
gravitational acceleration g may be neglected provided 



(2., 



Thus, for the air-water interface, gravity can neglected for wavelengths A = 27r/fc ^ 1.7 cm. 

In most experimental situations, w is a real, fixed frequency determined by the wavemaker and equation (|]^) defines 
the wavevector ko{Lo) as a function of u) for a particular a. When the viscosity of the fiuid is explicitly included in the 
linear analysis of surface fluctuations, the wavevectors have a positive imaginary componentlij (derived in Appendix 
A) corresponding to viscous damping of the surface fluctuations. 



B. Scattering by a circular surface domain 



We now consider the disturbance to the plane progressive wave caused by a circular domain on the interface. This 
domain changes the surface tension within it such that 



cr(r) aoQ{R -a)+ aiQ{-~R + a) 

where Q{x) is the Heaviside function which is unity for x > and is zero otherwise. 

The velocity potential in the lower half-space can written as an incident wave coming from X 
scattered velocity potential, '(/'(i?, 9, Z), which is assumed to be outgoing and to decay as i? — > oo: 



(b{R,9,Z,t) 



cosh ko{H+Z) 

( 

cosh koH 



,ikoX — iujt 



■ip{R,9,Z)e- 



(2.9) 
-oo, plus a 

(2.10) 



where fco is given by the dispersion relation (2.7) with a = ctq- Using the 9^—9 angular symmetry of the problem, 
the scattered potential may be expanded as 



tjj{R,9,Z) = cos m9 / dk%p^{k)J^(kR) 
m=o -^0 



cosh k{H+Z) 
cosh kH 



(2.11) 



where 9 is measured from the incident wave vector (see Figure la) and the Jm{kR) terms ensure that the potential 
remains bounded as i? ^ 0. Equation (2.11) is a complete expansion in Bessel functions and describes the analytic 
bmction ip{R, 9, Z) obeying the Laplace equation in the domain z < and satisfying the no-penetration condition at 
Z = —H. Later in the analysis, we will see that the choice of outgoing disturbances from the domain (a radiation 
condition) requires that the first kind Hankel functions, Hm\kR), represent the dominant contribution for large R. 
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The functions tpm{k) are determined by first substituting (2.10) into (2.6). Since the incident wave in (2.1C) 
automaticaUy satisfies the normal stress boundary condition at the interface, then outside the circular region (i? > a), 
where the surface tension is (Jq, we have 



['yoVidz + piu^MR. 0, Z)\z=o = (i? > a), 

while inside the circular domain we have 

[aiV\dz + pLo^]i){R,e,Z)\z=f) = kl{ai - cro)e''''''-^™"'' tanh fco^^ (0 < i? < a). 



We can further expand the right hand side of (2.13) using the identityc! 

oc 

Jo(fcoi?) + 2^iPJp(fcoi?) 



AkoRcosO 



■p=i 



(2.12) 
(2.13) 

(2.14) 



Substituting (2.14) and (2.11) into ( ^.13 ) and multiplying both sides of ( 2.13| ) and (2.12) by cos nO, then integrating 
6 from to 2tt, we obtain a set of dual integral equations for each n > : 



dk Jn{kR)il)n{k){puj'^ - fjofc^tanh kH) = 0, 

dk Jn{kR)ipn{k){puj'^ - (Tifc^tanh kH) = 

(2 - (5„,o)j"(o-i - cro)fco^n(fco^) tanhfco^^ 



(i? > a) 



{0<R<a) 



(2.15) 



Rescaling r — R/a,h — H/a,ka ~ q, and letting A — cjq/cti denote the surface tension ratio, we can rewrite this 
system of integral equations as 



Fn{q)Jn{qr) dq = 0, 
Fn{q)G{q)Jn{qr) dq = A^Jniqar), 



(r>l) 
(0 < r < 1) 



(2.16) 



where 



{q^ ta.nh qoh- q^tanh qh) . , . 
Fn[q) = :z i^n[q/a), 



(2.17) 



^_ q'^ tanhqh — qftanhqih 
q-^ tanh qh — qg tanh qgh ' 



and 



A„ = (2 - (5„,o)(l - ^)qo tanhqah. 
In ( |2l7| ), ( ^l8| ), and ( pl9| ), we have introduced the dimensionless wavevectors 



qj tanh qjh ■. 



(j-0,1). 



(2.18) 



(2.19) 



(2.20) 



We note that qi = kia denotes the dimensionless wavevector which would be measured if the entire interface were 
homogeneous with surface tension ai. 

Finally, as mentioned above, the wavevectors can be interpreted as having a small positive imaginary part which 
ensures that the potentials vanish as i? ^ oo. This behavior can be attributed to the effects of small viscous damping. 
The small damping may be explicitly accounted for by studying the viscous flow problem and taking the limit ^ 0+ . 
A dual integral equation formulation to the viscous flow problem for finite v is given for completeness in Appendix A. 
Here, we treat fully the inviscid problem, but simply allow the wavevectors to have an infinitesimally small positive 
imaginary part as is common in scattering calculations .Eil 
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C. Solution of the integral equations 



The flow problem contains three dimensionless parameters: qo^ the incident waveve ctor; h, the fluid depth; and A, 
the surface tension ratio ch aract erizin g the surface discontinuity. It is also clear from (2.19) that A = 1 gives A„ — 
for all n > so that using ( 2.16 ) and ( 2.17 ) 4'n = a nd hence there is no scattering. 

The solution of the dual integral equations (2.U) Sot the unknown functions Fn{q) can be readily determined 
numerically using a procedure introduced by Tranteij^ which we describe for completeness. Begin by noticing the 
form of the Weber-Schafhcitlin discontinuous integrals 



r(n + m+ l)r"(l -r2)'3-i 
q^^'^J2m+n+i3{l)Jn{Qr)dq^\ 2f> - {n + l)r {m + f3) 



TmiP + n,n + l;r'^) (0 < r < 1) 
(r > 1) 



(2.21) 



where !Fm{ + n , n + 1; r^) is the Jacobi polynomial of order m.c2l Based upon the analogy of the right-hand-side of 
( 2.21 ) and ( 2.16 ) for r > 1, it is natural to expand the unknown functions Fn{q) in terms of Bessel functions as 



F„(g)^qi-^^ A[:)j2™+„+/5(<?), 



(2.22) 



m— 



where /3 is a constant which may be chosen later to improve the ra tes o f convergence of the required numerical 
integrations. With this representation of Fn{q), the flrst equation in ( 2.16 ) is automatically satisfled. Substituting 
(2.22) into the second equation in (2.16), we are left with 



/ dqq'-^J2„^+n+0{q)Jn{qr)G{q) = A„J„(gor) 



< r < 1. 



(2.23) 



Multiplying both sides of (2.23) by r""'"^(l — r'^)^ + n,n + l; r^) and integrating r from to 1 gives for each n, 

where n = 0,l,2,3..., a system of equations 

^ / dqq'~^^J2^+^+p{l)J2l+nMq)G{q) = q^^ A^^J^i+n+M ? = 0, 1, ■ • ■ (2.24) 

Equation (2.23) represents an infinite matrix for the coefficients and can be written succinctly as 

00 

J2Ai"jL^^I{f3;qoA) = ql^'^A^J2i+n+p{Qo) {nj>0) (2.25) 

where 



m=0 



4"^(/3;9o,A) 



dx (x^ tanh qohx ~ A tanh qoh) 
^2/3-1 tanh (70/12; — tanh (70^) 
^2;-/m/^(go,^,A;x). 



J2m+n+f:i{qox)J2l+n+0{qox) 



(2.26) 



The infinite set of linear equations (2.25) determi nes th e coefficients of the expansion of the functions Fn{q) from 
which the Fourier coefficien ts ipn jk) are known by ( ^.17|) and hence the scattered wave ■(/'(?', ^, z) is determined. 



The integral in equation (2.2 



has a pole on the real axis at x = 1 (i.e. q — q^). We may treat this pole in a manner 
consistent with outgoing wavesHil by interpreting the pole as having a positive imaginary part, Im{qj} = e > and 
taking the limit e — > 0. This limit is consistent with viscous effects becoming small, — > 0^. In the limit e ^ 0, the 
integral can be evaluated by taking the path of integration through a small semicircle below the pole at x = 1 which 
yields 



- -(1 - A) 



sinh qf)h cosh q^h 
3 sinh go ^ cosh go ^1 + go 



J2m+n+l3{qo)J2l+'n+l3{qQ) 



+ I dx 
10 



I^^}{qo.h,K-x)+I^^}{qo,Kk-2-x) 



(2.27) 



dxI^iiqo,h,A;x) 
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where we have separated the integrable real part of the integral over from the contribution from the semicircle 
below the pole. The remaining integrals are evaluated numerically. 

We next solve a truncated system of linear equations ( 2.25 ) and so obtain the Am^'s which depend on go, h, and A. 
The scattered potential follows from 



^ cos n9 



cosh q{h + z) 



<7q tanh q^h — q^ tanh qh cosh qh 



J2ni+n+p{q)Jn{qr) 



(2.28) 



from which we find the interfacial height displacement of the scattered part of the wave 

q^~^ cos n9 



i i 
Vs{r,0) = -—'ip{r,9,z)\^=o = - 

U! OZ UJ 



qQ tanh qoh — q^ tanh qh 



J: 



f}{q)Jn[qr). 



(2.29) 



Equations (2.28) and (2.29) also have poles on the real axis at g = go, which we choose to integrate using a semicircle 
contour below the pole. This treatment fixes the scattered waves as outgoing. In particular, to evaluate the integral 
in ( ^.2^ ), we replace the J„(gr) term with ^{Hn\qr) + Hn\qr)). For r > 1, the integral over J2m+n+p{q)Hn\qr) 
is closed by the contour in the upper half-plane shown in Figure 2. The detour in quadrant one corresponds to a 
branch cut at g = a/ iwa^ jv and is derived in Appendix A. Similarly, the J2m+n+f3iq)Hn\qr) integral is closed in the 
lower half-plane. Only the Hn'' integral encloses the pole at g = go. Thus, only the outgoing waves proportional to 
H^\qQr) are important for large r. The final result can be expressed in the form 



dq 



cosh go (/i + z) 



gg tanh q^h — q^ tanh qh cosh q^h 



J2m+n+piq)Jniqr) 



cosh go /i cosh go (/i + 2 ) 
gQ~''^(3sinhgo/icoshgo/i -|- go^i) 



J2m+n+0{qo)HiPiqor) - r+ - r_ (r > 1) 



(2.30) 



The contributions to the integral from the integrations along the imaginary axis in Figure 2, F-t, can be shown to be 
negligible compared to Hn\qor) as r ^ oo. Asymptotically, 



HL'Hqor 



^iqQr — iTTn/2 — i'K/4: 



Trgor 



(r — > oo) 



which leads us to write the far-field potential in the form 

i/jir, 9, z; h) ~ J{9] A, h)—= — 

y/r cosh go n. 

where the scattering amplitude is obtained by combining ( ^.2^ ) and (2.30): 



(r oo) 



/(^;A,/i) 



/3+1 
% 



cosh go ft. 



3 sinh go ft cosh qoh + qoh 



E A^^'^ J2m+n+ 13 {qo) cos n9. 



n.m—0 



The total scattering section is denoted Zl and is given by 



I](go,A,ft)^/ d9\f{9)f = ^ 

% 



cosh go ft 



-1 2 



3 sinh goft cosh go ft + goft 



oo oo 



71—0 m— 



(2.31) 



(2.32) 



(2.33) 



(2.34) 



where denotes multiplying by 2 for the n = term. Equations (2.33) and ( ^.34 ) are the principle results of this 
section. 
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D. Numerical Evaluation of A 



(") 



The linear set of equations ( 2.25| ) was solved by numerically integrating the functions in ( 2.27 ) to obtain complex 



All numerical 



matrix coefficients and then, for each n, an AI x M linear system was solved to determine Ai-, 
integrations demanded relative error 10~*. The dependence of the Am'' 's on the matrix size was found to be negligible 
(< 10^^) provided M > 14. (3 was chosen to guarantee convergence of the integrals and we chose = 1 for 
our calculations as this appeared to give the fastest convergence for large values of m,L A numerical integration 
technique for infinite integrals involving products of Bessel functions, developed recently by Lucas,Lj was used for the 
integrations of the kernels in (2.27). 



E. Wave scattering on a heterogeneous free membrane 

Since we are not aware of an analogy to capillary wave scattering from a boundary condition inhomogeneity as 
studied above, it is useful to consider a related scattering problem with a similar geometry. In particular, we treat 
wave propagation on a stretched elastic membrane and calculate the scattering due to a circular region with tension 
Ti which is imbedded in a two-dimensional membrane with tension Tq. Vertical displacements of this membrane obey 
the usual wave equation. In Section IV the scattering from this model is compared with that derived in the previous 
section. 

The configuration described in this section can in principle be realized experimentally by embedding a thin flexible 
loop of thread in a soap film. Surfactant can then be added to the interior of the loop to decrease the surface tension 
within the circle. The thread now supplies the line tension 7, which acts in the tangent plane of the interface and 
satisfies 

l^{To~-Ti). (2.35) 

The effects of an underlying liquid substrate are not present in this example and the only dynamical variable of the 
membrane is the height fluctuation ?7(r, 9, t) which obeys a wave equation d^rj — c^'S/^rj with a dispersion relation 

2 2 

?' = ^ (i = 0,l), (2.36) 

where we have used the same rescaled dimensionless quantities as in the previous section (e. g., g = ka). The wave 
speeds for the different regions are Cj — \/Tj/ pm where pm denotes the areal mass density of the thin film. 

We now wish to consider scattering of an incident plane wave, e^wcosS-iwt^ this circular region of different 
tension. The membrane displacement, rye~*'^*, outside the circular region takes the form 

Vout{r,9) = e"i"''^°-<' +T^,{r,6) > 1) (2.37) 

where the outgoing scattered height can be expanded as 

7js{r, e)^J2 AnHl^\qor) cosnO. (r > 1) (2.38) 

n=0 

Similarly, the membrane displacement inside the circle can be expanded as 



,{r,e)^^BnJn{qir)cosne. (0 < r < 1) (2.39) 



At the domain boundary (r = 1), continuity of the sheet requires 

r]out{l,e)=ri„,{l,e), (2.40) 

which also implies dgi^ is continuous across r = 1. The second boundary condition requires careful consideration of 
the directionality of the line tension. Since both the interfacial tensions Tj and the line tension 7 act in the surface 
of the membrane, a balance of forces in the z direction requires continuity of slopes 

drr]out{l,e)=drllUl,e). (2.41) 
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Subs tituti ng the expressions for the height displacements inside and outside the loop into the two conditions (2.4C) 



and (|2.4lD yields 



A to X ^■« goJn(gi)[^+i(go) - J«-i(go)] - qiJniqo)[Jn+i{qi) - Jn-i(gi)] 



for n > 0. We note that the t ension ratio enters this result since qi/qo = ^Tq/Ti. At large distances, r oo, the 



membrane displacement (2.37) is 



r,„„t(r,0) ~ e"?«'- « + ^ (2.43) 



with 



FT' °° 

/(0;To/Ti) = J e-"/4 V(-z)"A„cosn0. (2.44) 

V '^qo 



n=0 

The scattering cross section for this process is given by 

4 9 °° 

E(go,To/Ti) = -l^oP + - V |A„p, (2.45) 
?o 90 

which will be compared in Section IV with the corresponding results obtained in IIC. 

III. SCATTERING WITH BENDING RIGIDITY DISCONTINUITIES 
A. Waves on a fluid coupled membrane with bending rigidity 

In this sect ion, we consider the case of a thin elastic plate overlying a fluid whose motion is assumed irrotational and 



satisfies (2.2). Imbedded in the plate is a circular region with different elastic constants. This circular domain thus 
acts to scatter an incoming transverse (bending) vibration. Physical manifestations of this problem are membranes 
with phase separated domains, with imbedded impurities or defects, or other films with a patch of different flexural 
rigidityt£l (see Figure lb). 

We assume that the on ly no rmal r estor ing force at the interface is an elastic resistance to bending which modifies 



the boundary conditions ( 2.12 ) and ( 2.13 ). Balancing the pressure, p = —pdt4>, the beiiding forces, and the inertia of 



the membrane with mass density pp ( |2.6D , the equation for the plate displacement C isE^ 

[Dv\ + ppd^^]C{R,e)^p. (3.1) 

Here D is the bending rigidity, which for an isotropic harmonic solid is typically defined aS 

- TTTTT^. (3-2) 



12(1-/1 



where d is the plate thickness, and Ej and pj are Young's moduli and Poisson's ratios, respectively. The elastic 
parameters will be assumed to have different, though uniform, values inside and outside the circular domain. In 
microscopic systems such as monolayers and bilayers where the "plate" thickness has a molecular scale, the ben ding 



rigidity is typically denoted k, but the physical interpretation of k in terms of material parameters as given by (3.2) 
is not valid. 

Taking a time derivative and using the kinematic constraint ( ^.3[ ), we obtain the dynamic boundary condition 

[DV\dz + {p + Ppdz)dl](t>{R, e, Z)\z=o = 0. (3.3) 

The solution of this problem is similar to that of the surface tension discontinuity. However, rather than considering 
the general case of arbitrary values of fluid density and membrane density, we simplify the problem and consider two 
limiting, important cases. The first case is when the normal forces from fluid pressure dominate the plate's own 
inertial response. The second situation, discussed in subsection B, neglects the effects of the underlying fluid. We 
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note an important difference between these two cases. In the formulation for the fluid coupled, massless plate, no 
further conditions on the behavior of the fields at r = l,z = can be imposed. Here, fluid velocities, surface 
displacements, and their derivatives are continuous. However, the physical model considered may require extra 
conditions at r = l,z = which cannot be satisfied by the velocity potential of the irrotational fluid. These extra 
boundary conditions are part of plate theory and are incorporated in the second case, the scattering of free plate 
vibrations. 



When the plate mass is neglected, the dispersion relation derived from (3.1) is 

Qj tanh Qj h 



(3.4) 



Equation ( p.4| ) is the massless plate limit, Pp/\ <C p, which, upon solving for A = 27r/fc using (3.4) for h = oo, is 
consistent with 



Pp < 



p'^D 



1/5 



(3.5) 



The detailed scattering calculation has the identical form to that of the surface tension discontinuity problem except 
that the functions F{q) and G{q) are changed and the coefhcients (2.26) are modified. It is straightforward to obtain 



■^w (90, Ad] x)= «7r(l - A^) 
f-i 



smhqQhcoshqoh 



dx 



Ssinhqo/icoshgo^ + 90^ 
4"? (go, h, Ad; x) + I^^} (go, h, A^; 2 - x) 



+l3{q{))J2l+n+p{qo) 

dxl';^i{qo,h,Ad;x) 



with 



r(")/„ u A . — ^1-2/3^^ tauhqphx - A^ tanhgpfe 



lZ'{qo,h,A;x)^x' 



x^ tanh Qohx — tanh qoh 



(3.6) 



(3.7) 



and Ad = Dq/Di denotes the ratio of out er to inner bend ing rigidities. The scattering amplitude / and total 
scattering cross section E t hen follow from ( ^.3S ) and (2.34) by solving for the necessary coefficients using (2.25), 
(with Ad replacing A), and (3.6). 



B. Bending waves of a free plate 



Here it is instructive to study scattering in the uncoupled plate limit. Provided appUp'^^ ^ p^/^I?^/^, the pressure 
from the underlying fluid can be neglected in the force balance (3.1) and it is no longer necessary to consider the 
underlying fluid (and so the solution for the velocity potential). Setting p = in (3T) yields a wave equation for C 
with the dispersion relation 



97 = 



Ppduj^a^ 



J = 0, 1 



(3.8) 



for waves outside (j = 0) and inside [j = \) the domain. We treat the dynamics Id a manner typically used in plate 
theory, where shear stresses and bending moments within the plate are balanced.EZI 

As in Section HE, the solution for ^ outside and inside the circular domain can be expressed as 



Co 



. iqQV cos6~iujt 



where 



and 



Cs(r,6) = [^n^i"(9or) + S„X„(gor) 

n=0 



(r > 1) 



cos n9 (r > 1), 



(3.9) 



(3.10) 



Qn{r,B) = ^ [CnJn{qir) + E^In{qir)]cosne (0 < r < 1). 



n=0 



(3.11) 
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To determine the coefBcients An, Bn, C„, and En, we apply four constraints at the boundary r = 1: continuity of 
the plate displacements, the pure bending approximation, balance of vertical shear forces, and balance of rotational 
moments. These quantities can be expressed in the following terms: — CmtT Pr — P^,M~^ = M+,, and — 
V^, where l3r,Mrr, and Vrz are the bending angle, the bending moment, and the Kirchhoff shcarj-stpess resultant 
respectively, and ± denote the functions evaluated as r — > 1 from outside and inside, rcspectivelyE3'c3 Expressing 
these four conditions in terms of the out-of-plane displacements, we conclude that 



D[r) 

or 



v^C + (MO-i)(r§ + iSl 



r2 dff^ 



(3.12) 



and 



D{r) 



dr 



la^ 

r 9612 



de^dr 



(3.13) 



are all continuous across the boundary of the circular inclusion. The latter two boundary conditions at the domain 
boundary are not enforced for the fluid coupled problem, as discussed in the next section. Substituting the expansions 
(3.9- 3.11 ) into these four boundary conditions yields a linear set of equations which determine the unknown coefficients. 
In particular, we find 



detN 
detM' 



(3.14) 



where the matrices M and N are given in Appendix B. Expressions for the scattering amplitude and total scattering 
cross section, which are functions of the An, are identical to those of the tension discontinuity problem and are given 



by equations (2.44) and (2.45) 



IV. DISCUSSION 



A. Results 



We first examine capillary wave scattering on a he terog eneous surface using equations ( 2.33D and (2.34), and 



compare with scattering in a free membrane, equation ( 2.45| ). Polar plots of the scattering amplitude, \ f{0, h = oo)|. 



for various incoming wavevectors go and A = 2 are shown in Figure 3. In the long wavelength limit, qo I, the 
domain is nearly a point scatterer, and the scattering is nearly isotropic as expected; even for go ~ Ij there is almost 
no variation with 6. Distinct lobe structures develop when the incoming wavelength is about three times the radius of 
the scattering domain. Figures 4a-d illustrate, respectively, the effect of subphase depth on scattering for wavevectors 
qo = 1.0,2.0,3.0, and 4.0, again for A = 2. In each case, depths h = 0.1,0.5, and h = oo are considered. For 
these wavevectors, decreasing h decreases the scattering amplitude for all 9, but has almost no effect on the angular 
variations. 

The total scattering cross sections as a function of the reduced incident wavevector for A = 2.0 and 10.0 are shown 
for various depths in Figures 5a,b respectively. The depth dependence of S is shown in Figure 6. Figures 5a-i) show 
oscillations in E as qo is varied. Unlike examples of sound wave scattering from an infinitely rigid obstacle,c3 where 
a monotonic variatioHL of S with go is expected, the calculations presented here exhibit a different and more complex 
dynamical behavior .Ej The curves for^S have oscillations and maxima analogous to Mie resonances in the scattering 
of radiation from a dielectric sphereE3 In Figure 6, a weak variation of S with h is observed which varies most 
significantly for intermediate go. We note that for shallow underlying fluids viscous dissipation is important when 

^ a?uih'^ , and is discussed further below. We have also compared our results with the^ymptotic expression in the 
small wavenumber limit (go ^ 1), S = ^(A — l)^g"^, developed by Chou and Nelson.tS The agreement is good as 
displayed in Figure 5c. 

We compare the above results with those of scattering in a free membrane (the effects of the underlying fluid 
are neglected, section HID), shown as the darker dashed curves in Figures 5a, b. Although these two problems are 
geometrically similar, they have different dispersion relations, so in order to provide a meaningful comparison, we 
study two cases with the ratio of wavevectors gi/go = 2^/^ and 10^/^ (A = cto/cti = 2, 10 but To/Ti = 2^/3, 10^/3). 



A plot of (2.45) for Tq/Ti = 2 would correspond to a scale change (compression) along the go axis in Figure 5. The 
dependence of E on go for scattering in the uncoupled membrane versus the capillary wave problem is qualitatively 
similar at the same ratio of inner to outer wavevectors. In fact, the scattering cross section for A = 2.0, /i = cxd in 
Figure 5a is nearly indistinguishable from that of the free membrane. The corresponding amplitudes \f{6,h)\ for 
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the uncoupled membrane, equation (2.44), are nearly identical to the polar plots shown in Figure 4 corresponding to 
A = 2.0. For A = 10.0, Figure 5b, the free membrane model is now clearly distinguishable from the capillary wave 
scattering case, even after the wavevectors are rescaled to account for the different dispersion relations between a free 
membrane and a capillary wave. 

To illustrate the effects of varying the ratios of the surface parameters, we have plotted in Figure 7 S((7oj A, ft, = oo) 
as a function of A^^ for go = 1-0, 2.0, and 3.0. The limit A^^ = (Ti/cto ^ oo is associated with an infinitely "hard" 
domain at the interface which resists increases in surface area, though we note that such monolayer systems are 
unstable and can exist only in the presence of an externally applied force. Increasing A from unity corresponds to a 
more flexible domain, which actively supports internal fluctuations, and leads to a more complicated scattering and 
far-field disturbance. 

We now present results for scattering in interfaces characterized by bending rigidity rather than surface tension. 
Scattering amplitudes for the fiuid-coupled membrane for qi/qo = 2^/'^ (corresponding to Dq/Di = 2^^'^) are shown 
in Figures 8a-d for = 1.0, 2.0, 3.0, and 4.0 and different fluid depths. Total scattering cross sections as a function 
of Qo are plotted in Figure 9. For comparison, the scattering amplitude and cross sections for bending waves of a free 
plate (no underlying fluid) are depicted by the dotted curves in both Figures 8 and 9. To maintain scale similarity 
in the internal to external wavelength ratio, qi/qo = 2^/-^, we choose Dq/Di = 2'*/'^. In addition, due to the elastic 
boundary conditions at r = 1 in the plate model, we need to specify the contrast in bothj-bending rigidity D and 
Poisson's ratio fi. The dotted curve in Figure 9 was calculated with fiQ = 0.25, /ii = — 0.25.Eil From other numerical 
simulations, we have found that the cross sections E are rather insensitive to the discontinuity in Poisson's ratios 
except for small variations in the flne structure. 

The discontinuities present in the higher derivative boundary conditions (3.12) and ( |3.13 ) for the free elastic plate 
make this problem mathematically different from that of the massless bending interface coupled to an underlying fluid, 
where the velocity potential is assumed harmonic. For the fluid coupled irrotational flow problem, only two boundary 
conditions at r = 1 can be satisfi ed. Solutio n of the full boundary value problem with four boundary conditions 
applied at the domain boundary, ( [3.12 ) a nd (|3.13| ), requ ires consideration of viscous effects in the underlying fluid. 
The discontinuities at r = 1 suggested by (|3.12D and ( 3.13 ) are smoothed out in the thin boundary layer with thickness 
y/v/uj adjacent to the upper surface. It appears that the differences in Figures 8 and 9 for wave scattering along 
the heterogeneous plate as compared with the comparable capillary wave scattering problem are the discontinuous 
boundary conditions at r = 1. In the surface tension problems solved in Section II, the membrane displacement and 
its first derivative are both continuous and so both (2.40) and (2.41) are satisfied. 

The oscillations in S(<7o) for surface tension scattering, bending rigidity scattering, and scattering in an uncoupled 
membrane or plate, (Fimures 5 and 9), are similar to those of the cross section of electromagnetic radiation scattering 
from a dielectric sphere.E2l In that case the radiation field penetrates a sphere with finite dielectric constant; internal 
reflections constructively or destructively interfere depending on the incoming wavelength. For the problem studied 
here, capillary wave scattering from an interfacial inhomogeneity, internal waves are supported, and interferences lead 
to the oscillations in T,{qo). 



B. Experimental Consequences 

Our analysis has been conflned to wave scattering from an isolated domain. However, under most experimental 
conditions the morphology of the monolayer film cannot be directly controlled. For example, when amphiphilic 
molecules are deposited at the liquid-air interface they form a film which is in a phase separated state consisting of. a 
collection of scatterers of different sizes. The density and sizes of these domains depend on surfactant concentration. H'EI 
Interfacial suefactant morphology thus affects capillary wave propagation with scattering as a source of enhanced wave 
attentuation.ErO |— ■ r—. 

Experiments on capillary wave dampingtiTtJ show at least one maximum in the attenuation coefficient as a function 
of total surface concentration of surfactant. Wave amplitudesi-^peasured by I77P, can be related to the scattering in 
the dilute scatterer limit (single scattering only) by Beer's lawE3 

|77(a;)|2 cx e-^^(*)^ (4.1) 

where c is the area concentration of domains of radius a. Since the calculated cross section E(go) contains peaks 
(Fig. 5a,b and Fig. 9) (tl)-|at least qu alita tively accounts for the maxima in the damping coefficients which are 
experimentally measured.E3nl3 Equation ( [4.l[ ) suggests that peaks may observed as go is varied. Indeed, the variation 
in damping as. domain size and incoming wavevector are tuned has been qualitatively observed in the experiment of 
Wang et. a/Jij who measured a sharp increase in damping as the incoming wavelength was decreased from A ^ a to 
A « a. 
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We conclude the discussion by examining the region of vahdity of the analytical results derived in Sections II and 
III. Throughout the analysis, effects of viscosity have been neglected. This approximation is valid provided the domain 
is larger than the viscous length, ^Jvju)^ i.e., the Stokes layer adjacent to the fluctuating boundary. Thus we require 




(4.3) 



When the dispersion relation (2.7) for infinite depth is used to eliminate to from (4.3), and a condition that the waves 



are not too heavily damped by viscosity is imposed, k~ -^lij/v^ we obtain validity criteria for qj 



2 \ 1/2 



> qj > 



2 \ 1/3 



(4.4) 



For an air-water interface, cr « 70 dynes/cm, v « 0.01 cm^/sec, and p «1 gm/cm^. The first relation in ( [1.4| ) provides 
an upper bound of kj — qj/a ^ cr/pv^ « 7 x 10^ cm"-'^ and a lower bound of kj = qj/a 3> 113 cm~^, 5.2 cm^^, and 
0.24 cm^^ for domain radii of a w 10/im, 100//m, and 1mm r espe ctively. 



For a rigid membrane described by the dispersion relation (3^), the conditions A;| uj/v and (O) lead to 



kja = qj 3> max 



av^p 



1/5 



aiP p 



9P_ 



1/4- 



(4.5) 



where the last terroJs a condition due to the neglect of gravity waves. For typical estimates of the bending rigidity 
in cell membraneSjO D ~ 10~^^ergs, the above condition (4.5) requires exceedingly small wavelengths, A <C Inm, to 
be satisfied. As an estim ate of the range of validity, consider a system with D = 0.1 ergs. For domain radii of 10/zm, 
lOO^m, and 1mm, (|4.5D requires kj » 63 cm^^, 10 cm^^, and 9.95 cm~^, respectively. 



V. SUMMARY 



In this paper, we have formulated an analytical method for studying surface wave scattering due to a circular 
inhomogeneity in the boundary conditions. The results of the calculations imply that scattering of capillary waves 
from a circular domain where there is a discontinuity in surface properties (surface tension or bending rigidity) is not 
very sensitive to the depth of the underlying liquid (see Fig. 7). The scattering patterns \ f{9\ h)\ have lobed structures 
indicative of angular variations in the far-field disturbance. In addition, an examination of the above results (Figures 
5a and 9) show that differences between scattering from a surface tension contrast and a bending rigidity contrast 
are small if the wavevector ratio, qo/qi, of the two problems is maintained constant. For a thin elastic plate, the 
scattering from a domain of different rigidity has some qualitative differences when compared with the scattering from 
the fluid coupled system (Figure 9); this difference is due to the additional boundary conditions applied at i? = a 
for the free plate problem whereas the neglect of viscosity requires the elimination of two boundary conditions in the 
irrotational flow problem. 

Of course, in real experiments of phase separated surfactant films, the domains are polydisperse in size, often 
interact, and may not be circular. To incorporate such effects, a statistical average over domain sizes and correlations 
in the scattered potential must be included. However, we expect the qualitative behavior of the simple models described 
here to be representative of the physical features of these systems. Furthermore, phase separating monolayers are rarely 
in equilibrium, the typical sizes and number density change after the ambient conditions are changed, and the time 
evolution of the film obeys complicated coarsening and topological relaxation mechanisms. It is not unreasonable that 
experimentally, long times need to elapse before reproducible damping results are attained. In fact, when scattering 
is important, capillary wave damping might be used as a diagnostic for probing the statistical averages and effective 
influences of phase separating morphologies. 

The calculations presented for scattering with irrotational flows can be extended to include systems of membranes 
completely immersed in a fluid. Common examples are cell membranes, emulsions, or bilayers. As long as the 
geometry is symmetric top and bottom, both side contribute equally to the dynamical pressure; thus, the normal 
stress boundary condition becomes (T(r)V^?y = 2pdt(i>. All the remaining steps follow with the substitution p 2p. 
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Results of scattering for such an immersed membrane are s imply those for scattering at an air-liquid interface with 



qo,qi rescaled by the modified (p —> 2p) dispersion relation (|2.20| ). 

We have only considered a circular domain which affects the normal stress balance at the interface. However, the 
mathematical methods utilized ia,this paper can be extended to study tangential stress variations due to viscoelastic 
properties intrinsic to the film.ll3 Furthermore, a similar approach may be applicable for studying the effects of 
nonuniform surface stresses on the propagation of surface acoustic waves. 
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APPENDIX A: VISCOUS EFFECTS 

In this appendix we reconsider capillary wave scattering from a circular surface tension discontinuity in the presence 
of substrate fiuid viscosity. We study solutions of the continuity equation V • u = and the time-dependent Stokes 
equation, dtv = —'S/{p/p) -\- v'S/'^v for small amplitude capillary waves in an infinite depth fluid. Here, all variables 
r, 2, and t are dimensional. As before, the system is driven at a fixed frequency u and all quantities are assumed to 
have an e~*'^* time dependence. The limit of small viscous effects corresponds to the analysis presented in Section II. 

When considering viscosity of the underlying liquid, tangential stresses at the interface must also be balanced. 
Accounting for the line tension 7 at the circular domain boundary and changes in surface tension cr, and neglecting 
any surface viscosity or elasticity of the surface film, the tangential stress balance has the form 

- -S{r -a)+ V_ia(r) + pv{d,v^ + V±_vA = 0. (Al) 
a 1^-" 



For a step change of surface tension treated here (equation 2.9), the discontinuity of V±cr{f) at r ~ a is balanced 



by the line tension term, -^Sir ~ a) and the surface tangential stress condition obeyed everywhere on the interface is 

{dzV^+W±v,)\^^^ = 0. (A2) 
The normal stress balance including viscous effects is 

a{f)Vlri = -p+ 2pvd:,vA^^ (A3) 



Taking a time derivative of (A3) and using (2.3), we have 

cr{r)V\vz + 2ipvujdzVz = iujp^^^^ (A4) 
To solve the unsteady Stokes equations, it is convenient to decompose the velocity fields a si 

z;=V0(r)-H(^^-^), (A5) 

where — 0, and u±_ and Uz satisfy the diffusion equation (recall all fields are proportional to e~*"*) 

— iuju±^ = i^V^iii — ibJUz = v\l'^Uz- (A-6) 

The pressure is p = iujpcj). 

We treat the incident (i) and scattered (s) waves separately by defining the incident fields, <j)^^\x, z) and (tij*'' , wi*^ ) , 
and scattered fields, ijj{r,9,z) and {u'i\ui^-'), such that 

</>(r, e, z) = (t)^'\x, z) + V'(r, e, z) 

= uj'^ (x, z) -I- (r, 6*, z) (A7) 
Uz = uk\x, z) + u-z\r, 0, z). 
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The incident fields (j){x,z) and {uf' ^u^f') determine the fluid velocity of a damped w ave on a uniform interface 
of surface tension ctq. These fields must satisfy the normal stress bou nda ry condition (A4) with a = and the 
two-dimensional divergence of the tangential stress boundary condition (A3), respectively, 



Co 



0. 



(A8) 



Consider an incident wave with a velocity potential of unit amplitude on the y-axis, x = z = Q. Substitution of 
0^')(a;,z) — e^Qo'^+Qaz ^^^^ ui\x,z) — ^e*'3o^+'oz where Iq — Qq — iuj/v, into equation (A8) yields two equations for 
the unknowns A and Qa{uj) which are satisfied when 



A- 



2vQl 



2vQq — iuj 



and 



(A9) 



(AlO) 



The condition (AlO) determines the dispersion relation for capillary waves^on a homogeneous surface with tension 
(Tq. In the ^ limit, the first correction to the pole at cro^Q 



uj'^a^ givesEj 



ko + ie = ko + i 



(All) 



which shows that incident waves are spatially damped by the factor gi-'^'^^/'i'^o)^ _ 

We now consider scattering from a domain of radius r = a and surface tension cri . The incident wave is the same as 
in the above nonscattering case, i.e., a wave that has unit velocity potential amplitude in the absence of the scatterer. 
The normal stress condition obeyed by the scattered fields is 



(r> 1) 
(0 < r < 1). 



(A12) 



Since the tangential stress condition holds for all r, we can use ( [A^ ) to eliminate the Uz^ terms from ( |A1^ ) hence 
solving for tp. 

Expansions for the fields '0 and (uj_,Mz) can be found which obey the Laplace and diffusion equations respectively. 
For each angular harmonic cosniO, ip{r, 9, z) can be written as a Hankel transform of il'rn{k)e'^^ as in equation (2.11); 



similarly wi*^ can be expanded in terms of Am{k)Jm{kr)e''^ . Furthermore, with u± — {ur, ug), it is convenient to work 



with uf'*'' ± which are expanded in terms of U^{k) Jm±i{kr)e''^ ; 

ip{r,9,z) = ^ coamO / dkt})m{k)J'm{kr)e^'' , 



(A13) 



?TI = 



u^^\r,e,z) = Yl '^o^™^ / * [t/+(fc)J™+i(fcr) + t/„(fc)J™_i(fcr)] 



(A14) 



and 



POO 

{r,e,z) = ^ sin 7^6* / dk [C/+(fc) J™+i(fcr) - [/,";;(&) J„j_i(fcr)] 

^ poo 

u[^\r,9,z) = ^ cosm6' / dk Am{k)Jjn{kr)e''' 

m=0 



(A15) 



(A16) 
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where I = -^/fc^ — iuj/i/ and we must now determine the four functions U^{k),Am{k), and V'm(fc)- Note the above 
integrands have a branch cut starting at fc = y^iuj/v. 

The r emain ing s teps foUow thos e de scribed earher though the algebraic details are more cumbersome. Substitution 
of ( A13 ) and ( |A16| ) into equation (A2), multiplying the integrand by cosnd and integrating 6 from to 2tt, we relate 
the two fields 



AJk) 



2vk^ 



2iyk^ - 



(A17) 



Applying the normal stress b alanc e on the interface in the interior of the circular scatterer yields the dual integral 
equations of the same form as (2.16) where now we find the dimensionless functions 



Fniq;Z) 



and 



G{q;Z) 



ql-q^ + AZ^q^m(q + m) ) i/'„(g/a) 



q^ - ql- AZ'^Kq^m{q + m) 
g-' — 5q — AZ'^ Kq'^ra{q + ra) ' 



(A18) 



(A19) 



corresponding to (2.17) and (2.18) respectively, where the strength of the viscous effects are measured by the Ohnesorge 
number 



(To a 



(A20) 



We have used the definition (2.2C) in the h 



limit for qj and defined the dimensionless quantity ni = la = 



q^ — qj"^ jZ. Also, the equivalent of the right side of the second equation in ( ^.16| ) is 
A„(Qo; Jn(Qor) = (2 - 5„,o)(l - A) 



1 + 2iQ2a2^/^- 



3/2 



(A21) 



Solving the the dual integral equations yields Fn{q: Z) from which the tpn{q/a) arc determined. When Z small, the 
capillary wave pole in G{q; Z) is g ~ Qq as given in ( |AllD . This positive imaginary shift when ly approaches zero fixes 
the integration path below the pole at q w fco in ( |2.26| ). In the — > (Z — > 0) limit, this problem reduces to that of 
scattering of surface waves on an irrotational fluid of infinite depth. 



The scattered velocity fields can be extracted by using the decompositions of the scattered fields, (A13-A16) to 
balance the f and 9 components of the tangential stress equation (A2). The remaining expansion coefficients U^{k) 
are solved in terms of the ipn{k): 



Ut{k)^± 



vkH 



2l/fc2 



ILO 



-^n{k). 



Adding the contributions from 7/'(r, 0, z), the scattered velocities can be formally written as 



oo ^ 

(r, 9, z) — ^ cos 1719 I 

m=0 



dk 



2vkl 



2vk^ 



dr 



(A22) 



(A23) 



4'^(r,6l,z) 



m=0 



sin m9 



dk 



2vkl 



2vk^ 



— Jm{kr)%prn{k) 
r 



(A24) 



{r,9,z) — cosm9 / 



kdk 



kz 



2vk^ 



2vk^ 



Jyn{kr)il),n{k). 



(A25) 



To obtain the total velocity field of the scattering process, the incident velocities, (corresponding to 0*^*-' and 



[uf' ^ui"^) in (|A5|)) must be added 
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% ~ 6" 

Iq 2vQq — iuj 



(A26) 



— n^p'^Qax+Qoz _ '^^Qo iQox+loz 



{k27) 



where Qq is given by the roots of ( |AlC| ).lli 



APPENDIX B: BOUNDARY CONDITIONS FOR BENDING WAVES 



This Appendix tabulates intermediate results for the scattering produced by a circular inclusion with bending 
rigidity Di imbedded in a membrane with bending rigidity Dq. The four boundary conditions formulated in Section 
III B lead to four linear equations for the coefhcients, Xi = (A,i, _B„, C„, En) which may be compactly written 



B, 



where 



B = i^{2 - 5„,o) 



Jniqo) 



and M is the 4x4 matrix given by 



(Bl) 



(B2) 



(B3) 



M2^ = f-i/(i)'(go),-K(go), -J'M, -IM) 
\ 90 go / 



(B4) 



M31 
M32 

M34 

M41 
M42 

M43 
M44 



(9o - fJ'Qn'^)Kn{qo) + qo{^J.o - l)^n(9o) 

^ [{ql + (Ml - 2)n2)J„(gi) - 91(^1 - l)4(gi)] 

-s^ [(9? - fiin')in{qi) + - mAqi)] 

-q,{ql - Mon2)i?i^^'(9o) - n2(Mo - l)i?i'^(<7o) 
'Zobo + (/"o - 2)n2]x;(go) - n2(Mo - l)if„(go) 
^ [(7i(g2 _ ^W).r,^{q^) + (mi - l)n2j„(gi)] 

[91(9? + (Ml - 2)^2)/; (91) - (mi - l)n2/„(9i)] 



(B5) 



J'n{x), etc. denote the derivatives of Bessel functions. 

The matrix N is identical to M except that the first column, Mu, is replaced by the column vector B. 
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FIG. 1. (a) A membrane with a circular domain of surface tension ai imbedded in an interface of tension ctq . (b) A plate 
of thickness d and flexural rigidity Do with a circular inclusion of rigidity Di. Ripples in the surface represent the membrane 
fluctuations due to an incident plane wave and the scattered wave. 



FIG. 2. The contour used to evaluate Lj^l and 4>{r,6,z). An infinitesimally small positive viscosity and imposition of an 
outgoing solution requires that either the pole requires a small positive imaginary part, or that we take the contour through a 
semicircle below the pole. 

FIG. 3. Scattering amplitudes \f{0)\ for an infinite depth system with A = ao/ai = 2. The same plots for the uncoupled 
membrane, equation (2.50), at To/Ti = 2^/^ are nearly identical. 



FIG. 4. The scattering amplitudes \ f{9;h)\ for various depths h = oo, 0.5, 0.1 at go = 1,2,3,4 for A = 2. A "+" identifies 
the origin. 

FIG. 5. (a) Total scattering cross section E(q(),A = 2,h) for the surface tension discontinuity problem. Four depths are 
considered: h — oo, 0.2, 0.1, and 0.01. (b) Total scattering cross section E(go, A = 10, h) for h = oo, 1, and 0.1. (c) A comparison 
of the scattering cross section calculated numerically and asymptotically in the small qo limit for A = 1.1. In both (a) and (b) 
the dark dashed curve is E in the free membrane case. Section HE. 



FIG. 6. Depth dependence of E(go, A = 2, h) for scattering from a circular surface tension variation at incoming wavevectors, 
qo = 1.0, 2.0, 3.0, 4.0, 5.0, 8.0, 9.0, and 10.0 

FIG. 7. A dependence of E(go, A,h = oo) on surface tension ratios. For A^^ > 1, mechanical stability requires an externally 
imposed negative line tension. The limit A~^ — > oo corresponds to an unstretchable domain. 

FIG. 8. The scattering amplitudes for an interface with bending rigidity; h = oo, 0.5, 0.1 and A^ = 2*^^. The dotted curves 
are scattering from a free elastic plate. In this case, the Poisson's ratios are fj,o = 0.25, /Ui = —0.25. 

FIG. 9. Total cross sections of bending wave scattering for fluid depths h = oo, 0.2, 0.1, 0.01. The dotted line is for a free 
plate with fj,o = 0.25, /xi = —0.25. 
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